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^\ Abstract 
O ■ 

' The perturbations in the electron number density during recombination contributes to the cosmic 

C/3 ! 

microwave background bispectrum through second order terms. Perturbations in the electron 

, density can be a factor of ^ 5 larger than the baryon density fluctuations on large scales as 
> ■ 

' shown in the calculations by Novosyadlyj. This raises the possibility that the contribution to the 

' bispectrum arising from perturbations in the optical depth may be non-negligible. We calculate 

, this bispectrum and find it to peak for squeezed triangles and of peak amplitude of the order 

OO ' of primordial non-Gaussianity of local type with f^i ^ 0.05 ~ —1 depending on the ^-modes 

o ■ 

' being considered. This is because the shape of the bispectrum is different from the primordial one 

^ . although it peaks for squeezed configurations, similar to the local type primordial non-Gaussianity. 
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I. INTRODUCTION 



First order perturbation theory has been of sufficient accuracy for analysis of the Cosmic 
Microwave Background (CMB) observations so far. However future CMB experiments will 
have high enough precision that second order effects would need to be taken into account 
for theory to have similar accuracy. The second order contributions will in particular be 
important for the higher order statistics like the three point correlation or the bispectrum. 
Second order effects in CMB have been studied previously 
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13| . Bartolo et al. have derived the Boltzmann equations at second order and also 



the analytic solutions for the CMB transfer function at second order with some simplifying 

nn n 

assumptions [IJ, |15|, see also [16 1. 

All numerical and analytic calculations at second order so far have ignored the contribu- 
tion arising from the perturbations in the electron number density, 5e- These contributions 
are expected to be small compared to other second order terms since 6e multiplies the col- 
lision term which has contributions from the difference of first order radiation and electron 
dipoles , radiation quadrupole and higher order moments of the radiation transfer function. 
These terms are small during recombination compared to the monopole terms. Recombina- 
tion however depends on matter and radiation densities and perturbations in the electron 
number density can be quite different from the perturbations in the matter and radiation 
densities. This was calculated by Novosyadlyj who showed that this is indeed the case 
and perturbations in the electron number density can be ~ 5 times the baryon number 
density perturbations. 

We calculate the CMB bispectrum on all scales arising due to the perturbations in elec- 
tron number density and compare it with the bispectrum expected from a primordial non- 
Gaussianity of the local type. This is a full numerical calculation without any other ap- 



proximation except that we only consider terms involving 6e- Although this 



bispectrum 



turns out to be below the detection levels of future experiments like Planck 18|, there are 



some important general implications which are discussed in the conclusions section. We use 
the following cosmological parameters for our calculations (values at redshift z = unless 
specified): baryon density = 0.0418, cold dark matter density Qc = 0.19647, cosmo- 
logical constant Q/^ = 0.76173, number of massless neutrinos N^, = 3.04, Hubble constant 
Hq = 73, CMB temperature Tcmb = 2.725, primordial Hehum fraction i/He = 0.24, redshift 
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of reionization Zri = 10, primordial gravitational potential power spectrum P{k) = 2ti'^ /k^ 



II. INHOMOGENEOUS RECOMBINATION 



|l9l| by Novosyadlyj [l7|, which is a modification of the 
20| to calculate the perturbations in the electron number 



We use the code DRECFAS^ 
recombination code RECFAST 
density 5^ = {rif. — ne)/ne during recombination, rif. is the local electron number density and 
Tie is the mean electron density. Perturbations in baryon (5^) and photon density {6-y) result 
in perturbations in the electron density with an amplitude that is amplified or suppressed 
depending on which terms in the evolution equations prevail. Specifically photo-ionization 
prevails on superhorizon scales resulting in (5e ~ 5 x 5^ during recombination. We refer the 



reader to 



17( 1 for further details. 



We will not consider the full second order Boltzmann equation [14)] but only the terms 
involving the perturbed electron density. This is given by: 



dr 



-k5p 



9(1) 



+ n.Vb - ip2(H.n)n« 



where r is conformal time and tq its value today. O = AT/T = O^^^ + O^^^ + 
higher order terms is the fractional perturbation of CMB temperature, superscripts indi- 
cate the order of perturbation while subscripts denote the multipole moment. All other 
perturbations are of first order and we will omit the superscript for them. Vector quantities 
are in bold face and their magnitudes in normal face with " denoting unit vectors. We have 
omitted the factor of 1/2 usually multiplied with the second order term [l^ for convenience. 
n is the unit vector along the line of sight, k = dn/dr = —rieaTCi is the mean differential 
optical depth due to Compton scattering, ctt is the Thomson scattering cross section and 
a is the scale factor. We take the electron velocity to be equal to the baryon velocity V5. 
V2 is the Legendre polynomial of order 2 and 11*^^) = 02^^ + BpQ + is the polarization 
term, subscript P denoting the polarization field [2l|. We must caution that this partial 
equation is gauge dependent because 6e depends on the gauge. We will be using conformal 
Newtonian gauge for 6e. The combinations of terms multiplying 6e is gauge invariant. All 
perturbed quantities are functions of r and coordinates on spatial slice x. G is in addition 
a function of line of sight angle h. 



Following standard procedure 
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22l | , we take the Fourier transform of Equation [T] and 
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integrate formally along the line of sight. 



^0 



TO 



(2^ 



Mk-k',T) 



&^o\k',T) -e^^\k',h,T) + n.k'Vb{k',T) - ^V2{k'.h)U^^\k\ 



T 



.(2) 



where g{T) = — /t(r)e '^^^^ is the visibility function and k(t) = JJ" dT'ne{T')aTa{T'). Also 
Vij{k', r) = k'Vb{k', r). We now take the spherical harmonic transform of Equation [2]to get 

(2) 

the multipole moments, Q^rn- 



Q^'\k,n,ro)YUn)dn 



This integral can be performed after decomposing B*^^-* into multipole moments, 
6^^)(fc',n, r) = J2i"{—iY" i"^^" + l)Vi"{n.k')&^^) {k' , t) and using relations between expo- 
nential, spherical harmonics, spherical Bessel functions and Legendre polynomials 23|]. Note 
that which is the dominant term in the multipole expansion of B*-^-*, will cancel out. 

We will see later that the dipole term partially cancels the effect of (H), the Vishniac term. 
So only i > 2 modes in B'-^^ which are expected to be small compared to monopole, will 
contribute to the bispectrum. The result is: 

e?l{k,ro) = rdrgir) [ ^,6^{k-k\T) 

JO J [271) 



e'm'e"^Om" 



1(4» 



i'm'm" 
2 



\ 47r(2^+ 1) 



^ - /^jj-(2£ -\- 1) ^^^'V'0^£'m'£''m''if'[^('^ '^o)]Ye'm' i^)^e"m" i^')'^e" 



(2i" + 1)3 
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{2i' + 1)5 



(3) 



d^k' 



drgir) I --^5e(fc - k', T)S'"'{k, k', k', r) 

J (ZTTj 

CeJ^'e"m" Clebsch- Gordon coefficients, je are spherical Bessel functions. The sums are 
over all allowed values of im with the exceptions explicitly specified. The last line defines 
the function S*^™". Its arguments are written so that we can keep track of the part, k', that 
statistical variables like temperature anisotropy depend on from the part that deterministic 
functions depend on, k', separately. 
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III. BISPECTRUM 



We can now use Equation [3] to calculate the bispectrum. This is defined as: 

B^mrmlm, = i^fXi ^o)^^^ ^o)(^f3m3 ^o)) + 2 pemiutatious, (4) 

where aimix, Tq) are the coefficients in the spherical harmonic expansion of the corresponding 
temperature anisotropy. () denotes the ensemble average. At second order they are just the 
Fourier transform of 0^^(fc, tq) while at first order they can be computed from <d^l\k, tq). 

We can now calculate the first term of the bispectrum in Equation HI 
(1,1,2) = {afXA^,ro)a^lJm2i^:ro)afJ^^{x,To)) 

~ ^^""^ J f27r)3f27r)3f27r)3^ ^ ^hn.Af^^)^e,mAf^2) 



(27r)3 (27r)3 {2tt) 

W) 



fj drg{r) J ^,{Se{ks - k', T)S'^"^^{k„ k', k', r)Q'i^{k^, ro)eg)(fc2, Tq)) 



(5) 



We can write each term in the ensemble average as a transfer function times initial gravita- 
tional potential perturbation. Thus, 

Se{ks-k',r) = <^,{ks-k')6,{\ks-k'\,r) 

eg)(fci,ro) = Uki)Q^l^ih,ro) 

S^^"'^{ks,k',k',T) = $i(fc')5^""'(fc3,fc',A;',r) 

(<l>,(fci)<l>,(fc2)) = {2nf6%k^ + k2)P{k^) 

We are using same symbols for statistical variables and their deterministic transfer function 
counterparts, with arguments determining which one we mean. ^3 jg f^]^Q three dimensional 
Dirac delta distribution and P{k) = 27r^/A;3 is the initial power spectrum. Since we assume 
the initial perturbation to be Gaussian, the 4-point ensemble average can be decomposed 
into 2-point ensemble averages. 

($,(fc3 - k')Uk')^,{k^)Uk2)) = ($.(fc3 - k')^,{k'))){Uk^)<!>,{k2)) 



+ ($i(fc3 - fcO$^(fcl)))($.(fe')$.(fc2)) 
+ ($,(fc3-fc')*.(fc2)))($.(fe')<^'.(fcl)) 

+ i2TTf6%ks + fci - k')P{h)5%k2 + k')P{k2) 

+ (27rf 53(^3 + - fcOP(A;2)53(fci + fc')P(A;i) (6) 

First term in Equation [6] contributes only for fcs = 0, it is a product of monopole and 
power spectrum and is unobservable. The other two terms are identical with fci,fc2 terms 
interchanged. So we need consider only one of these. Denoting the two terms by superscript 
(1, 2) and (2, 1) we can write the first term of the bispectrum as: 

(1,1,2) = (l,l,2)(i'2) + (l,l,2)(2.i), 

{27i)^6\ks + fci - k')P{h)5\k2 + k')P{k2) 

Pik,)P{k2)5,ik,,T)S''^'{ks, -k2M.T)Qf^{k,,T,)Q^l^{k2,To)5\k^ + fc2 + fcs) 

(7) 

In the last step we have used one of the Dirac delta distributions to integrate over k' . To 
proceed further we use the representation of Dirac delta distribution as Fourier transform 
of unity and the expansion of exponential function in spherical harmonics. 

5\k^ + k2 + ks) = J ^6*^*=^+'"^+'"^^" 

e^C--) = A7rY.^'jiikr)Y;^{k)YUr) (8) 



Using Equations [8] in [7] we can perform all angular integrals and all radial integrals except 
two which involve transfer functions of perturbations and the line of sight integral. The 
integrals involving spherical harmonics result in Wigner 3jm symbols which can then be 



23|. The result after performing these 



summed using, for example, formulas tabulated in 
integrals is: 




(9) 



TT J 



/'+^2+«^_l)€2(2f' + i)(24 + 1) 



^ f'2 ^2 f" 



(f^2,r)je' [k2(To - r)] 







+m{k2,T)j',^[k2{To - t)] + -n^'\k2,T) {SJ^Mto - t)] + j,, - t)]} 



- / kldk2P{k2)efJ{k2,To) 
TT J 

_ ^ /'+^2+^2(-l)^2(2r + l)(24 + l) 
j1[A:2(t-o -r)] 



^ ^'2 ^2 f" 







+ [eb{k2,r)-e^{k2,r)] 



+ -n(i)(fc2,r) {SjIMto - r)] + J,, - r)]} 



(10) 



In the last step we have defined iVb = Ob/k and 6^ = 3kQi and evaluated the sum over £'2 
explicitly for 0i. It can be seen from this expression that the effect of Vishniac term 61, is 
partly cancelled out by 6j. In this form the gauge invariance of Bqq is also apparent. In 
arriving at these expressions we have also used the identity je{—x) — {—lYji{x). The prime 
on the Bessel functions denotes the derivative with respect to the argument. 

We can now write down the final expression for the angular averaged bispectrum defined 

by: 



/ 



B 



m\m2raz 

V 47r \^ y L 



(11) 
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IV. PRIMORDIAL NON-GAUSSIANITY OF LOCAL TYPE 



We will compare our results with the bispectrum from primordial non-Gaussianity of 



local type. This is given by 
2 



24 



25j: 



B: 



prim 



(2£i + l)(2£2 + l)(2£3 + l) [ h 4 ^3 





47r 



Piir) 



2 




(ir(ro - 



edkP{k)Qeik,To)ji[k{ro-r)], 
k^dkfNLQi{k, To)ji[k{To - r)], 



(12) 



where f^i is the non-Gaussianity parameter defined by the following form for the primordial 
potential, $i(a;) = $L(a;) + /atl ($|(a;) - ($|(a3))) with $L(a;) Gaussian. Note that the 
expression for Pe is similar to B^q and Bqq, the difference being the additional modulation 
by the terms at recombination in the later case. As we will see later, is similar in shape 
to the visibility function g{T) but peaks at an earlier time. All plots and results are for 

fNL = 1. 



V. NUMERICAL CALCULATION AND RESULTS 



We calculate Se in conformal Newtonian gauge using DRECFAST 171]. All other first 
order terms are calculated using CMBFAST 2^. In particular Q£ii{k, r) is given by the line 
of sight integral: 

e^»(A;,r) = e"(") f dr' S^^\k,r')3e>[k{T - t')] (13) 
Jo 

Here S^^\k,T') is the usual first order source term. Since we are evaluating the transfer 
function at r < tq, we g et an extra factor of e'^^'^\ otherwise this is same as the standard 



line of sight formula 22] . Qe" becomes smaller with increasing i" and we cut off the sum in 
Equation [TO] at £" = 30. This is accurate for r < lOOOMpc which is sufficient for the present 
calculation since the visibility g{T) is non- negligible only for 240Mpc ^ r ^ SOOMpc (Figure 
[2]). Wigner 3jm symbols are calculated using the code by Gordon and Schulten 2^ which 
is publicly available at SLATEC common mathematical library 271 ]. 

Figure [1] shows a comparison of Pi^r) and B^q{t). The modulation by Se results in 
shifting the peak to later times. Also visibility g{T) can be compared to primordial term 
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FIG. 1: /5£(t) and B^q^t) is shown as a function of r for several values of i. 
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FIG. 2: (tq — r)^a£(r) for several values of I and the visibility function g^r) as a function of 
conformal time r. (tq — r)^a£(T) peaks earlier than g{T). 
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FIG. 3: Bqq{t) is shown for several values of £. Also shown are contributions from the polarization 
term H, slip term 9b — 9g and from all the other terms ^i>2 '^^P ■ 
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FIG. 4: 0.1 X i{l + l)/3^(n),0.01 x i{£ + l)S^e(n), i{i + 1)BI^q{t^) and contributions to it from 
polarization, slip and rest of the terms is shown as a function of multipole moments L Some of the 
functions have been scaled as specified above. 
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FIG. 5: Absolute value of B^l^^^' labeled "Primordial" and S^^^^^a labeled "Recombination" for 
^3 = 10. Z axis is on linear scale while color plot shows the same on log scale. 

«£. They are similar in magnitude but have a different shape (Figure [2]). Bqq is however 
much smaller in magnitude than the other terms, Bse and jSe, as can be seen from Figures [3] 
and m at low i but become comparable at high i. This results in a much smaller bispectrum 
from recombination at low i compared to the primordial one. Figures El O [7] and [8] show 
the absolute value of the bispectrum from the primordial non-Gaussianity with /tvl = 1 and 
that due to inhomogeneous recombination for £3 = 10, 200, 1000, 2000 as a function of £1, £2- 
Z-axis is on linear scale while the color map is on log scale. They are almost identical at 
the peaks but differ considerably away from the peaks which occur when either ii or £2 is 
equal to £3 and the other is small, a signature of the local nature of the non-Gaussianity. 
At low i the amplitude 5^i^2^3 jg much smaller than -Bp^^^ but they become comparable at 
high i. Their signs are however different and this will become apparent when we estimate 
the confusion in /tvl due to B^^^'^^^ . 



To estimate the confusion to the estimate of fjqL we follow |28| and define the statistic 



Q prim 

^rec — 2^ n n n 

- hL E 77f^ (14) 

The result of solving Equation [14] for /at/, is shown in Figure [9] as a function of imax^ where 
Iraax IS the maximum value of I included in the sum in Equation [TH As expected from the 
examination of bispectra, j^i^ is small and positive at low i^ax but ~ — 1 at high imax- 
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l3=200 




FIG. 6: Absolute value of B^l^^^' labeled "Primordial" and S^^^^^a labeled "Recombination" for 
^3 = 200. Z axis is on linear scale while color plot shows the same on log scale. 
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FIG. 7: Absolute value of SpJ.-^^' labeled "Primordial" and S^i^^^a labeled "Recombination" for 
I3 = 1000. Z axis is on linear scale while color plot shows the same on log scale. 



VI. CONCLUSIONS 



We have calculated the CMB bispectrum due to inhomogeneous recombination. This was 
expected to be small because the combination of terms multiplying 6e is small. However 



calculations by Novosyadlyj 17| showed that 6e could be large and this suggested that 
the CMB bispectrum could be non-negligible. Although it turns out to be small it is still 
larger than what one might have expected from making an estimate based on tight coupling 
or instantaneous recombination approximation jisj and ignoring the perturbations due to 
inhomogeneous recombination. This is especially evident at high imax- Also the bispectrum 
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FIG. 8: Absolute value of B^l^^^' labeled "Primordial" and S^^^^^a labeled "Recombination" for 
^3 = 2000. Z axis is on linear scale while color plot shows the same on log scale. 
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FIG. 9: Comparison of primordial bispectrum from local type non-Gaussianity with bispectrum due 
to inhomogeneous recombination in terms of parameter /tvl as a function of imax, the maximum 
£ mode considered. 



from recombination looks remarkably like the local type primordial bispectrum, which is not 
entirely unexpected since both arise due to product of two first order terms. Since the other 
second order terms in the Boltzmann equation [l4| are expected to be larger than the ones 
we considered, our calculation motivates a full second order numerical ca 
terms in order to assess their effect on future experiments such as Planck 



culation of these 
and the level 



to which they cause confusion when probing for primordial non-Gaussianity. 
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